Demônio de Creutz Implemente o “demônio de Creutz” para o modelo de Ising, na rede quadrada. Compute as médias da magnetização e da susceptibilidade para vários valores de energia e diferentes configurações iniciais. Compare os resultados com o Q2R, colocando os resultados no mesmo gráfico. Compare o tempo de relaxação da magnetização neste método com o Q2R. Encontre a relação da temperatura com a energia: a) plotando o histograma da energia do demônio $E_d$ e b) através da relação $kT/J = 4/\ln (1+4J/\bar{E})$

Percolação na rede quadrada